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We compute the rates of real and virtual photon (dilepton) emission from dense QCD matter 
in the color-flavor locked (CFL) phase, focusing on results at moderate densities (3-5 times the 
nuclear saturation density) and temperatures T ~ 80 MeV. We pursue two approaches to evaluate 
the electromagnetic (e.m.) response of the CFL ground state: (i) a direct evaluation of the photon 
self energy using quark particle/-hole degrees of freedom, and (ii) a Hidden Local Symmetry (HLS) 
framework based on generalized mesonic excitations where the p meson is introduced as a gauge 
boson of a local SU(3) color-flavor group. The p coupling to generalized two-pion states induces a 
finite width and allows to address the issue of vector meson dominance (VMD) in the CFL phase. 
We compare the calculated emissivities (dilepton rates) to those arising from standard hadronic 
approaches including in-medium effects. For rather large superconducting gaps (several tens of 
MeV at moderate densities), as suggested by both perturbative and nonperturbative estimates, the 
dilepton rates from CFL quark matter turn out to be very similar to those obtained in hadronic 
many-body calculations, especially for invariant masses above M ~ 0.3 GeV. A similar observation 
holds for (real) photon production. 

PACS: 12.38.Bx, 13.40.Hq, 12.40.Vv 


I. INTRODUCTION 


Recent extensive studies of QCD at finite quark chemical potential (^q) have established that the ground state at 
large densities and zero temperature is a color superconductor with large pairing gaps on the order of 100 MeV (see, 
e.g, ref. 0] and references therein). Specifically, for 3 flavors of massless quarks, general arguments based on free 
energy minimization and maximal symmetry favor the formation of a condensate that effectively locks rotations in 
color and flavor space Among the intriguing implications of this color-flavor locking (CFL) are the modifications 
of photon propagation in such matter, e.g., true electromagnetism being a combined rotation of the standard U(l)^ 
and U{1)y color hypercharge. The color-flavor locked medium is thus transparent to this rotated photon [Q. More 
recent studies have addressed the question of reflection and refraction at the interface of normal (pairing-free) and CFL 
media |^. Given a finite (small) temperature, a natural question would be to ask for the blackbody radiation of the 
system, i.e., how much a block of CFL matter shines via emission of real photons or dileptons (virtual photons). For 
photons, the answer is pertinent to their propagation in the cores of neutron stars, which may be sufficiently dense for 
formation of color superconducting quark matter. A study of the temperature dependence of the photon emissivities 
and mean free paths in the CFL phase can determine how transparent this phase is at various temperatures. 


Apart from consequences for specific physical systems like neutron stars, one ought to address the issue of photon 
and dilepton emission at varying temperature and baryon density as a general means to investigate different regions 
of the QCD phase diagram, especially in colder and denser regions. The latter may be accessible to the HADES 
experiment at GSI . Small final state interactions make photon and dilepton probes important long lived signatures 
of the state of strongly interacting matter. Corresponding production rates have been extensively studied for both hot 
hadronic matter incorporating medium effects ill (see r ef, for a recent review) as well as for the Quark-Gluon 
Plasma using perturbative [p^[-p0[ and nonperturbative [^-23 approaches. It has been pointed out that perturbative 
qq annihilation rates in a QGP at comparable temperatures and densities are suggestive of a ’quark-hadron duality’, 
i.e., the dilepton signature for an interacting confined phase with broken chiral symmetry resembles that of a QGP 
with chiral symmetry restored. This has been found for both the low-mass {p-/u}-) region |24| and the intermediate- 
mass region [ p5| , p6| (between the (j) and J/5'). Gould such a dual scenario hold for the dense phase of condensed 
diquarks at finite temperature as well? In this paper, we will address this question by making perturbative estimates 
of the photon and dilepton emissivities from dense superconducting matter with 3 flavors of massless quarks (the GFL 
phase). While such perturbative calculations are strictly valid only at parametrically large chemical potential, it is 
intriguing to extrapolate down to densities of a few times that of normal nuclear matter (ng = 0.16 fm“^). In this 
way, our approach attacks the finite density problem from the side of large baryon density and serves to complement 
studies of emission rates from the low density (hadronic) side, extended to include medium effects. Significantly, our 
calculation at finite density allows for finite temperature effects that are possibly as large as the gaps in the CFL 
phase, conjectured to be several tens of MeV at /i, ~ 300-500 MeV. 
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Another issue we have sought to address in our analysis of the dilepton emission rates is the viability of vector 
meson dominance (VMD) at high density. Along the lines of ref. we have included excitations describing the 

rho meson (p) in the CFL phase via an effective theory based on hidden local symmetry. Based on weak coupling 
analyses to leading logarithm accuracy, effective Lagrangian approaches to QCD in the CFL phase have proven to be 
useful with applications to color-flavor anomalies, hidden local symmetry and meson properties |^,^,^. Here, we 
use the effective Lagrangian to obtain the coupling of the p meson to the rotated photon in the CFL phase and model 
its width in the dense phase via a Breit-Wigner resonance through its coupling to two generalized pions. In this way, 
we recover VMD and show that the p meson makes an important contribution to the dilepton emission rate, much 
like in low density hadronic frameworks. 


The organization of this paper is as follows. In section we compute the photon polarization tensor from the quark 
loop in the CFL phase (independently confirming the result first obtained in ref. |]^ ), and present numerical estimates 
for corres pond ing dilepton and photon rates at finite temperature and for baryon densities down to n_B = (3-5) no. 

we derive the effective Lagrangian for the CFL phase which incorporates hidden local symmetry and 


In section III 


show how the generalized p meson arises as a dynamical gauge boson of this symmetry. In particular, we address the 
fate of vector meson dominance in this model, thereby determining the coupling of the generalized vector mode to the 
photon. Again, we present the associated contribution to the dilepton production rate. In section IV we confront the 
combined emissivities from CFL matter with prior estimates from in-medium hadronic as well as perturbative QGP 
calculations. We summarize and conclude in section]^. 


II. CFL PHASE AND E.M. CURRENT CORRELATOR I: WEAK COUPLING 


A. One-Loop Photon Self Energy 


The photon polarization tensor in color-flavor locked quark matter to one-loop order was first derived by Litim 
and Manuel ||33t] . For the sake of completeness, and as an independent confirmation of their result, we here provide 
a calculation using a different quark basis than in ref. |^. In what follows we adopt the notation of ref. Q]. We 
assume all quarks to be massless. The self energy receives contributions from the diagonal (dressed quasiparticles) 
and off diagonal (diquark) components of the Nambu-Gorkov propagator which reads ||^ 


cf^^_^G+(A) S-(A) 

>^1^1 - S+{K) G- {k) 

where (G±)-i = GJ - E±, with $+ - (V'c^), 

denotes the charge conjugate partner of tp. The off-diagonal components of the matrix in eq.([l]) are ^ = —"-rg 


( 1 ) 


^(AT) = Hgjo. ipc 


{iptpc) and (GJ) \ . 

he matrix in eq.(||) are = — G(j^<i)^G^. 
With the convenient choice of SU{3) generators T® = Xs/^ = {\/3/2)Q, {Q = diag{2/3, —1/3, —1/3}) and T® = 
diagjO, 1/2, —1/2}, the mixing of the photon with the eighth gluon can be described by the linear combinations 


Ai_i = Af^cosOcFL — G^^sinOcFL 
G® = A^sinOcFL + G^ cos 9c FL 


( 2 ) 

( 3 ) 


where tan = 2el\/Zg, and e, g denote the gauge couplings of A^, G^ respectively. The GFL quarks couple to 
the in-medium photon A^ with strength e = e cos 9c fl via the vertex 




iX8y^Sab-6^^{Xs)ab r 


( 4 ) 


since Q = Q^l + l^Y where Y represents the hypercharge operator. In eq. (^), {i,j} and {a, b} denote flavor- and 
color-indices, respectively. The self-energy then takes the form 


= -^Trfc,s,c./,7VG 


T^S{k+^)T^S{k-^) 


( 5 ) 


where the traces are to be carried out over internal 4-momentum (^, spin (s), color (c), flavor (/) and Nambu-Gorkov 
{NG) indices. The NG trace is easily performed in view of eq. (P, whereas the color-flavor locked structure of the 
condensate is conveniently unraveled using the projectors ]3q], 
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( 6 ) 


1 


1 


CL ^ fisi cL ^ - 5isi), cL ^ ^{SabS^^ + sisi) - cL . 

These projectors may be identified with the singlet (Pi) and octet (Pg) ones according to 

Pi=C\ P8 = C^ + C\ 


( 7 ) 


Using the following trace relations, 


Ti-cj 

Ti-cj 


QPiQPi 

QPsQPs 


= 0 = TycJ 

= 4 = -Tr.j 


QPiQPs 


= Tr 


C,/ 


Q(C'"-C'^)Q^Pi 


Q{C^ - CLQ {C^ - cL 


the self energy of the CFL photon can be expressed in compact form as 


Uf^^{K) = 2e^ 


GZ(K)+^Z(K) 


( 8 ) 


( 9 ) 


with 


GW) = Y. - ^)) (10) 

e=± ^ 

2^8 W = E - f)) . (11) 

e=± ^ ■' 


This result for explicitly agrees with the one quoted in ref. cf. eq.(27) therein. The subscript 88 indicates 
that only the octet projection contributes. The sum over e = ± encompasses the contributions of particles (+) and 
antiparticles (-). As the antigap dependent pieces arise only at order 0{p,~L we neglect pairing effects for 
antiquarks in our calculation of emission rates. Moreover, the antigap is gauge dependent even for forward scattering 
of antiquarks We keep, however, the contribution from free antiquarks. This, in particular, has the advantage of 
naturally recovering the dilepton rate from the perturbative qq annihilation process for vanishing quark gaps, which 
is our baseline for assessing the effects of pairing. The Dirac spin traces may be performed explicitly using equations 
(45) and (95) of ref. p^ . The 4-momentum trace over k involves a 3-momentum integration as well as a summation 
over Matsubara frequencies in the energy component fco within the imaginary time formalism to incorporate finite 
temperature effects. This allows to obtain explicit expressions for the various components of n^'^(Ar) (see ref. ||^ for 
details), which in the following section will be employed to compute photon and dilepton emission rates. 


B. Photon and Dilepton Rates 


The medium modification of the photon fieldQ in the CFL phase is contained in the components of the polarization 
tensor, which in turn can be related to e.m. emissivities |^,^. The thermal emission rate of real photons reads 

^ ~(2^ = q) f^iqo;T) , ( 12 ) 

while that for dileptons (virtual photons) is given by 

W ' /hw;r), (13) 

where f^{qo',T) = l/[exp((j'o/T) — 1] is the thermal Bose factor, a = I'^t: with e the usual electromagnetic charge, 

and = q^ — q^ the invariant mass of the lepton pair. The factor of a implies that we choose the projection 


'^We refer here to the (massless) rotated photon in the CFL phase which plays the role of true electromagnetism. 
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along the ordinary photon which couples to the dileptons. The imaginary part of the retarded photon self energy, 
= Imn^+2Imn^, is the weighted sum of its transverse and longitudnal projections (at the photon point, i.e., 
for vanishing invariant mass = Qo ~ vanishes). It is obtained from by performing the appropriate 

analytic continuation to arbitrary complex values of the external frequency iK^ = go- The one-loop result for positive 
go > 0 then reads 


= r’'‘ 

ei,e2=± 


2 V eie2 


ki • k2 - 1 


S{qo - ei - £2) (1 - IVi - N2) + [( 5 (go - ei + £2) - S{qo + ei - £2)](IVi - N2) 


^{qo ~ £1 ~ £2) (1 — iVi — N2) ( 1 — ^ ^ ^ -I- [^(go — £1 + £2) — iJ(go + £1 ~ £2)] (-^1 — N2) 1 -I- ^ 


£i£2 


£i£2 


, ( 14 ) 


where the quasiparticle energies are denoted by ei = \J , with = Ciki — fj,q and ki = \k\, k 2 = \q — k\ the 
quasiparticle momenta. Ni = f^{ei\T) = l/[exp(ei/T) -|- 1] are the usual Fermi distribution functions which become 
©-functions in the T=0-limit. The summation over all combinations of ei = ± and 62 = ± generates a total of 4 
contributions embraced in curly brackets. However, since we approximate antiparticles as free (ungapped), the 
term is zero except for the ei = 62 = -I-. The physical interpretation of the various terms is encoded in the delta 
functions which determine the on-shell kinematics of the process (in eq. (p^, we have omitted terms involving delta 
functions of the form (5(go -|- £1 -I- £ 2 ), as they do not contribute to the production of on-shell photons or lepton pairs). 
Let us focus on the contribution of the particles to ImH. At zero temperature, Pauli Blocking enforces the distribution 
functions to vanish identically, i.e., Ni = N 2 = 0. Thus the only contribution to the imaginary part arises from the 
(5(go —£1 —£ 2 )-terms. For vanishing gap (free quark gas), energy-momentum conservation in the particle-hole excitation 
implies = — 2 fcifc 2 (l + cos0) {9 = Z(fci, ^ 2 )), i-e., the response is strictly spacelike which corresponds to standard 
Landau damping (this obviously remains true for T > 0). For a finite gap, the response sets in at go > 2A, the break¬ 
up of a Cooper pair acting as the threshold for excitations, as expected. More importantly, the response extends 
into the timelike region as well, which again follows from energy-momentum conservation: choosing, e.g., cos 6 = — 1 
(which is the optimal configuration for a timelike response), one finds = 2 (A^ -|- eie 2 — C 1 C 2 ) to be positive, thus 
allowing for photon and dilepton emission. These features are confirmed in Fig. |^, where Im n(go,g, A;r = 0) is 
plotted versus energy at fixed 3-momentum g = 60 MeV and for various values of the gap parameter. 

At finite temperature, Ai, A 2 7 ^ 0 and the delta functions with energy arguments go T £1 4= £2 also contribute. The 



FIG. 1. ImH at temperature T = 0 MeV, three momentum g = 60 MeV and for A = 0 MeV (solid), A = 10 MeV (dotted), 
A = 30 MeV (dashed) and A = 50 MeV (dash-dotted). In each case, the response has a threshold at go = 2A. 

resulting response function as shown in Fig. || indicates a significant relocation of strength from the spacelike to the 
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FIG. 2. Imll at finite temperature T = 50 MeV, three momentum g = 60 MeV and for A = 0 MeV (solid), A = 10 MeV 
(dotted), A = 30 MeV (dashed) and A = 50 MeV (dash-dotted). The relocation of strength from the spacelike to the timelike 
region is evident. This signals the possibility of dilepton emission. 


timelike region. Such a feature was also found in ref. where neutrino scattering in the color superconductor has 
been studied. The delta functions with opposing signs on the quasiparticle energies are restricted to spacelike PP, and 
are relevant for scattering of CFL photons off quarks (Landau damping in the CFL phase). In our context of photon 
and dilepton emission, where > 0, these terms are of no direct concern. However, we may note that the shifting 
of response strength from the spacelike to the timelike region implies the suppression of scattering as compared to 
the free Fermi gas [ ^ . 

To display the associated dilepton emission rates, it is common to focus on the three-momentum integrated quantity 
which produces the invariant-mass spectrum 


dRee _ f d Q dRee 

dM'^ J 2qo d'^q 

The pertinent numerical results for the CFL phase based on the one-loop photon self energy, eq. (0), are displayed 
in Fig. 1^, where only the particle contributions are retained to highlight the impact of Cooper pairing on the emission 
rate. The curves have been evaluated for a quark chemical potential of = 350 MeV (corresponding to a baryon 
density ub — 5.4no including massless strange quarks) assuming a temperature dependence of the gap as inferred 
from BCS theory, 


A = 




(16) 


The rather large variation in the considered gap values (10-70 MeV) is motivated by the sharp curvature of the above 
relation near p. For definiteness, using Tc ~ 0.56Ao with Aq — 150 MeV at /ig ~ 350 MeV (as suggested by both 
perturbative and nonperturbative g] estimates), the gap increases from a few to several tens of MeV over a small 
range of temperature around T = 80 MeV. The CFL rates are compared to free quark-antiquark annihilation at order 
Og in a baryon-rich QGP first computed in ref. 0: 


dRqq^ee _ P T g • Ti V (a^- + exp[-(go -k ^lq)/T]){x+ + exp[-fj,q/T]) 

d‘^q dTT^ g ql^ds^'^ "^{ 3 :++exp[-{qo + fj.q)/T]){x-+ exp[-fj.q/T]) ' 
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with a = e^/47r and x± = exp[—(go ± g)/2T].^ As characteristic for soft emission rates generated from in-medium 



FIG. 3. Dilepton rates at finite gap and finite temperature (with particles only). The rate vanishes as the gap goes to zero. 
The solid line is the 0{a^) estimate of the free qq annihilation rate at the same temperature and density. 

effects, the gap-induced processes exhibit a much steeper slope than the perturbative qq annihilation, exceeding the 
latter towards small masses. It is straightforward to incorporate the leading (free) antiparticle contribution within the 
CFL one-loop results by retaining the antiparticle projectors in the diagonal parts of the Nambu-Gorkov propagator. 



FIG. 4. Dilepton emission rate including the antiparticle contribution.(For A = 0 MeV, the perturbative qq annihilation rate 
is recovered.) 

Upon doing so we recover the limit of perturbative qq annihilation as the gap goes to zero, cf. Fig.^, where also 
the additional processes due to finite pairing gaps have been included. Whereas the gap acts as a threshold energy 
(go > 2A) in the response function, eq. (0). this is not so for the invariant-mass spectrum, as there is no constraint on 


^One should note that perturbative corrections arising from new processes allowed in a QGP at finite baryon density can 
give substantial enhancement over the lowest-order (©(og)) result of eq. ( |l7| ) in the low-mass region (M < 1 GeV) region [^ . 
Qualitatively, such corrections have a similar shape to the ones induced by finite gaps. 
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the three momentum transfer q. Nevertheless, the energy threshold leaves its traces also in the invariant mass spectra. 
To illustrate this point, we have magnified in Fig. ^ the region of very low masses M < 0.2 GeV. An inflection point 
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FIG. 5. Dilepton emission rate in the very low mass region. A smooth turning over seems to occur at M = 2A. 

of the rate is discernible at M ~ 2A, below which a finite 3-momentum transfer is required, entailing a suppression of 
dilepton production. In fact, the rate with the smaller pairing gap even exceeds the one with the more tightly bound 
Cooper pairs. 

The photon emission rate, ecu (p^, from the CFL phase at the same density and temperature (fig = 350 MeV, 
T = 80 MeV) is plotted in Fig. |, for different choices of the gap value.In the forthcoming sections, we will compare 
these rates to that from a hot hadron gas with in-medium effects included. 



FIG. 6. Direct photon emission rates from gapped quarks in the CFL phase at pairing gaps of 10 and 70 MeV (full and 
dotted lines, respectively). 


III. CFL PHASE AND E.M. CURRENT CORRELATOR II: HIDDEN LOCAL SYMMETRY 
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A. Generalized Mesons and Effective Lagrangian 


Beyond the 1-loop result from theprevious section, the existence of mesonic bound states in the vector and axial 
vector channel in high density QCD |^| is expected to further affect the dilepton emissivity from dense low temperature 
matter. The knowledge of the coupling of these modes to the electromagnetic current affords us the opportunity to 
complement studies of dilepton emission from hot and dense hadronic environments with “hadronic-type” rates from 
the high density side. Starting from a suggestion made in ref. (^, we will here adopt the Hidden Local Symmetry 
(HLS) approach to construct the coupling of these so-called ’generalized’ vector mesons in the CFL phase to the 
photon field. As an effective low energy theory, the HLS framework determines the couplings between generalized 
vector and pseudoscalar mesons, the latter being the Goldstone modes of chiral symmetry breaking in the CFL 
phase. This feature also provides the answer to the question whether dilepton production through the generalized 
vector mesons bears any overlap with the quark one-loop description of the previous section. It turns out that the 
quark-based (weak coupling) analysis, carried out to leading logarithm accuracy, does not support a coupling to the 
generalized pseudoscalar states (though it may be finite at next-to-leading order) and moreover, the coupling of the 
vector excitation to the physical vector current is zero Therefore, the decay of the generalized vector mesons 
(or, equivalently, the annihilation of generalized pseudoscalar mesons) adds to the dilepton rate independent of the 
one-loop result of the previous section. The HLS framework allows us to determine the coupling of the generalized 
vector mesons to the electromagnetic current, and hence the contribution to the dilepton rate. 

The idea of interpreting low-lying vector (and axial vector) mesons in free space as dynamical gauge bosons of a 
hidden local symmetry originated in refs. [|l3[-[l6[ . 

A particular advantage of this framework becomes apparent in the conventional [SU{2)l C S'C/(2)/{]giobai nonlinear 
sigma model where the electromagnetic interaction can be implemented as a global gauge symmetry, coupling to global 
isospin charge, distinct from the (hidden) local SU{2)v charge to which the rho meson (p) couples. This separation of 
the photon source from the p allows for a natural explanation of several phenomenological facts (such as p dominance 
of the vector coupling to pions, KSRF relations and universality of the p coupling) at the cost of fixing one parameter 
(the universal vector coupling constant). 

In weak coupling, the (composite) generalized vector mesons in the CFL phase are distinct from the screened 
and Higgsed gluons |^. Could they be the realization of a hidden local symmetry, besides the explicit local color 
symmetry? If the Cooper pairs have finite size, there is, strictly speaking, no hidden local invariance in the order 
parameter apart from explicit color symmetry. However, we will work in the zero size approximation, in which case the 
condensate does indeed exhibit a hidden local symmetry corresponding to the unbroken color-flavor SU(3) group^. 
Just like the p meson in case of hidden local symmetry at zero density |45|, the vector mesons of the CFL phase 
are composite and heavy. In this respect, our hidden gauge approach is different from ref. wherein the local 
SU{3)c was the hidden (color) gauge. In the following we construct an effective Lagrangian for the Nambu-Goldstone 
bosons arising from the spontaneous symmetry breaking, with the vector mesons introduced as gauge particles of the 
unbroken local symmetry. As in the vacuum case, they are rendered physical degrees of freedom by postulating a 
kinetic term (assumed to be generated at the quantum level), which does not affect the low energy dynamics of the 
theory. 


To start the derivation of the HLS Lagrangian in the CFL ground state, we recall that, formally, a nonlinear sigma 
model based on the coset manifold G/H is gauge equivalent to another model with Ggiobai C iJiocai symmetry and 
that the gauge bosons corresponding to the hidden local symmetry in iJiocai are composite fields. For the CFL phase, 
G = [S'17(3)c C SU{3)l C *S'C/(3)i{]giobai and H = [5'C/(3)c'-i-L+it]iocai- In terms of (right-handed) quark fields, the 
explicit form of the order parameter reads 


{<IrA) = , (18) 

where {a, 6}, {j,j} and {a,/9} denote color, flavor and spinor-indices, respectively (C: charge conjugation matrix, 
Pr\ right-handed chiral projector). An equivalent form applies to the left-handed fields. Following refs. 
this symmetry-breaking pattern implies the emergence of 8 Goldstone bosons corresponding to iS'17(3)c+A-valued 
condensate rotations of the form e"*'’'®® ^ with associated group generators {A = 1,...,8). In what follows 


^ A similar caveat holds for effective field descriptions of hadrons in the nonperturbative vacuum; it may be overcome by higher 
order terms in the chiral perturbation series. 
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we will assume the CFL ground state to be invariant under the parity transformation. We introduce unitary helds 
corresponding to the individual left- and right-handed rotations, characterized by expectation values 


{^Lax) = —{^Rai) = 5ai , ( 19 ) 

which implies ki = —K 2 in eq. (|T^), and hxes the overall U{1) phases in 

^l,r = Ul.r{x) ( 20 ) 

at, e.g., {(j)L) = 0 and {(j)R) = 7r/2 0 - The generalized left- and right-handed SU{3) helds can be parameterized by 

[/^(a;) (21) 

Ur{x) = , (22) 

where 'Kl,r = r^^^ the SU{3) generators are normalized as Tr(T"^T^) = 6^^ /2. At this point we have 

doubled the number of colored (generalized) Goldstone bosons represented by SU{3) phases of Ur and Ur. The latter 
transform under G as Ur ^ 9 tURg\ and Ur ^ respectively. The following further decomposition, 

Ur{x) =iJ(.x)iR{x) 

Ur{x) =^J{x)^r{x) , 

then makes explicit the additional invariance under the hidden local symmetry h G [S'17(3)c'+L-i-ij]iocai: 

^lix) ^ ^c(x)h^(x) (25) 

^lm(x) h(x)^R^R(x). (26) 

We will hx the local SU(3)c+r+r gauge by choosing (x) = ^r,r(x) = ^(x). It follows that ttr = ttr, so that we 
recover the expected number of Goldstone bosons (eight). This simply means that color locks to the havor diagonal 
so that left- and right-handed ’pions’ are the same. In this way, we implement the description of (generalized) pions 
as phases of the (colored) diquark condensate, and avoid an explicit appearance of the gluon helds in our effective 
theory^. With this realization, one introduces a vector gauge held Vp, that transforms inhomogeneously under the 
hidden local symmetry parametrized by h(x). Since the gauge helds may be viewed as connections in the coset 
manifold, we can dehne the gauge covariant derivative as 


(23) 

(24) 


• ( 27 ) 

The effective Lagrangian for the GFL pions is now assembled from invariants under SU{3)c G SU{ 3 )r 0 SU{ 3 )r G 
[5'C/(3)c-i-L+fl]iocai G (parity): 

^eff = - e ?^)' - ^Tr(A?t • C - C^)" ( 28 ) 

_ ar^TviDo^^ . ^ + Do^ ■ - 2iVof - as^Tr(AC^ • ^ + - 2iV,f + ..., 

where the breaking of Lorentz invariance of the hnite density ground state is reflected by individual temporal and 
spatial constants Ft and F’s, respectivelyj^ The hrst two terms of eq. ( p8|) reproduce the kinetic term for the pions 
while the other two (proportional to ur and as) vanish identically utilizing the equation of motion for the held, 

y; = -^Tr{T“(a^^^t + a^^t.^)} . ( 29 ) 


will not include in our low energy description the breaking of U{1)b —> Z 2 associated with baryon superfluidity (lY-nieson) 
or the anomalous t/yi(l) current ( 77 ' meson), as they are not relevant to the problem at hand. 

®In other works, it was the field E = U^Ur which mapped the meson octet, with gluons acquiring a mass via the gauge 
fixing of color as the hidden local symmetry. This integrating out of gluons led to an identification with the effective chiral 
Lagrangian at zero density. 

®One has Fs/Ft = = 1/y^ where Vtt is the velocity of the Goldstone modes [^. To leading logarithm accuracy, the pion 

decay constant Ft cs p^,po|-^. 
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As in the zero density case, the vector field propagates by postulating a kinetic term assumed to arise from interactions 
on the microscopic QCD level. This does not affect the limit of low energy dynamics that the effective Lagrangian 
is aiming to describe. The ellipsis in eq. (p8| ) denotes higher order and possible mass terms of heavier excitations 
consistent with the CFL symmetries, which we ignore. With the addition of the gauge kinetic term, and expanding 
up to 0(j}^) (p: momentum of the pion field), we have 

^eff = + Y(i9j7r)^ + ^Vo • x doiv + ^Vi ■ n x d^Tv 

+ l{aTF^)Vo^ + + ■ • • . ( 30 ) 

where we define the standard field strength tensor = 9^14 — d,^V^ + x Vj^. Upon rescaling the vector fields, 
Vo ^ grVo, Vi gsVi, we identify 

gVoTT-n- = -^ ar gr , = ar i ( 31 ) 

= -^ as gs , mg = as gl Fg . ( 32 ) 

We note that, apart from Lorentz symmetry breaking, the effective Lagrangian, eq. (^), is formally identical to that 
obtained at zero density with hidden local symmetry for the chiral SU{2)v subgroup. For the numerical calculations 
presented below, the value of gx will be taken to be that of gpTrTr in vacuurqj. We will show later 4-dimensional 
transversality of the photon self energy requires gs = gx- Both the mass of the vector field and its coupling to 
two generalized pions depends on the parameters ax and as- A suggestive choice will emerge when we introduce 
electromagnetism into our effective theory. This is our next step, so as to assess the coupling of the generalized 
(neutral) vector mesons to the (modified) photons with application to e.m. spectra. 

The true electromagnetic gauge boson in the CFL medium is the modified photon (A^), with a modified coupling e 
to charged degrees of freedom. is the gauge field associated with local U{1)q phase transformations on the Ul.r 
fields. Noting that this is equivalent to a simultaneous transformation ^ ^^^-igYsme^ 

that Ux,R transform as 


Ulm ^ , (33) 

where we have used Q = —Y for multiplication from the left. Thus, the covariant gauging with electromagnetism can 
be achieved via 


^rO=5^e/+*ee/ig (34) 

(35) 

(the superscript * denotes complex conjugation). Thus the effective Lagrangian, eq. (|^), with electromagnetism 
incorporated becomes 


nem 

^eff 




em^2 _ ^ J An 


y/Sgi 


V^gs 


-( 


ax — 2 


)eAo/AB87r^9o7r^ - ( 


as — 2 


2 3g^ 

)eAifAB8T^^diTT^ 


1 Ae^rrig r 2 

- x^Ai 

2 3g| 


(36) 


It obviously exhibits generalized p °-7 mixing (U® is assigned the generalized neutral p meson). Vector meson domi¬ 
nance (VMD), i.e., the absence of a direct TTTry vertex, is recovered for ax = as = 2, in analogy to free space, where it 
successfully reproduces the empirical p dominance of the electromagnetic pion form factor. Such a parameter choice 
thus seems natural on account of the formal resemblance of the CFL and vacuum case at the level of the bare effective 


^Here, boldface notation indicates 8-component vector in flavor-space with the cross- (x) and dot-products to be read as 
(^ X V)a = fABcV^V^ and V ■ V = with /abc the antisymmetric structnre constant of the SU{3) group. 

^Another estimate of gx may be obtained by a study of the 4-pion vertex by integrating out the p fields in the infinite mass 
limit. 
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Lagrangian. On the other hand, it has recently been argued based on a one-loop renormalization group equation 
(RGE) analysis, that VMD (for Nf=3) is only accidentally realized in the chirally broken phase at zero density, and 
that it does not hold towards the phase boundary of chiral restoration. The CFL phase, however, is a priori of a rather 
different nature (i.e., its symmetries and physical excitations), so that VMD cannot be precluded by the perturbative 
arguments invoked in ref. 

With VMD enforced, the mass eigenstates of the photon and the vector meson fields follow from eq. via 
diagonalization, 


Bfj. — {g^ + 

K = (9^ + (38) 

to| = 0 , = 0^2(5^-k ^e^) , (39) 

where g = gxiF = Ft for /i = 0 and g = gsiF = Fs for ^ = i. We see that p-y mixing spontaneously breaks 
[SU{3)]hidden 0 V(1)q down to U{l)em- 

We have thus shown that the HLS framework, which has proven successful in the chirally broken vacuum state, 
can be implemented into the (chirally broken, iVy=3) high density color superconducting ground state as well. The 
(generalized) p meson emerges as a resonance of generalized two-pion states, with its mass and width determined by 
the effective Lagrangian eq. (^^. We are now in position to evaluate its contribution to the dilepton emission rate in 
dense CFL matter. 


B. Dilepton and Photon Rates 


Adopting the VMD hypothesis for the coupling of (real and virtual) CFL photons to the generalized isovector 
mesonic current (as motivated in the previous section), the electromagnetic correlation function can be saturated by 
the p° meson. The pertinent dilepton rate, eq. (p^, then takes the form 


dRee 

d^q 


37r3M2 


af^[qo-,T) 


{mrY 


Poo l-m.D'^'^{M,q) 


ImD*^(M,g) 

9s 


(40) 


with a = /Fk, fP the thermal Bose-Einstein occupation factor, and the in-medium propagator of the generalized 
p meson, which is summed over its longitudinal and transverse polarizations to obtain the rate. The form of this 
propagator follows from inverting the equations of motion for the spatial and temporal components of the field 


Di^‘' = 




py 


q^q’' 


-^S M2 - ml 


% 

S 


2 _q±_ 


"It 90 


m'gq^ 


(41) 


Here, Pl and Pt are the standard longitudinal and transverse projectors for photons in matter |^]. We note that 
for mr = nis = riip, this propagator reduces to that for the massive spin-1 field in the covariant (Landau) gauge, 


= 


1 

M2 — m2 




(42) 


Though we have made a particular gauge fixing, the gauge parameter would appear only in the q^q'^ term of eq. (|^), 
and therefore not contribute to the rate by transversality of the photon. The rate is a gauge invariant (observable) 
quantity, driven by the imaginary parts of the longitudinal and transverse modes, i.e., spectral functions. The 
transverse propagator has a pole at M = ms, while the longitudinal propagator admits two positive definite solutions. 
These are given by 


n ml 
- — 


l±Jl-4q' 


. m^ — m' 


2 1 


(43) 


As mr > ms (cf. footnote 6), there exists a limiting 3-momentum given by gmax = mlg/{2^m\ — m?g) beyond 
which the longitudinal mode is damped. The lower mode corresponding to the choice of negative sign in the above 
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equation vanishes at q = 0 where it has zero strength, and only the upper mode persists, becoming degenerate with 
the transverse part, as it must in the rest frame. 

We may view the longitudinal mode as having an energy-momentum dependent self-energy by decomposing the 
inverse propagator as 

(44) 

Si(M, q) = - {ml - ml) (45) 

We impose retarded boundary conditions by continuing the 1/M^ factor as qo ^ qo — ie. To illustrate the above 
features, we insert a small (finite) imaginary part into the self energy and display real and imaginary part of the 
longitudinal propagator in Figs. and respectively, for fixed values of the masses {ms = 500 and tot = V^ms)- 



q„(MeV) 


FIG. 7. The real part of the longitudinal p propagator as a function of frequency for different 3-momenta q at fixed masses 
ms and mr- The upper and lower modes are degenerate at g = 5max and separate with decreasing q. 


For q approaching (jmax (— 176.5 MeV with our choice of masses), the two modes merge while for smaller q they are 
distinct. For q — 0, only the mode degenerate with the transverse part survives. The retarded boundary condition 
renders the spectral function, 


Al{co) = -2l-mDL{uj) 


(46) 


a positive definite quantity (which directly figures into the e.m. rates), even though the lower (’-’) mode has a negative 
discontinuity. We have checked that the spectral sum rule for the two longitudinal modes, 



du! 

2tt 


ll!Al{uj) = 1 


(47) 


holds for independent values of tot and ms- 

We can go beyond the mean field results, eq. (0), by including loop corrections from the coupling of the p to 2 
pions. In that case, the imaginary part of the p self-energy diagram (decay widtW are ascribed to the ptttt interaction, 
which we evaluate according to the ptttt matrix element determined by eq. (^). The self energy is generated to 
one-loop order from the squared matrix element integrated over available phase space and summed over final states, 

- ImE^^(go, Cl) = \Y. J ’ (48) 
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q„(MeV) 


FIG. 8. The imaginary part of the longitudinal p propagator for fixed ms, rriT and various values for the 3-momentum q. 
The upper mode is dominant except when q is very close to (?max = 176.5 MeV. 


with the reduced matrix element T = igpTrTre^(2q + k)^, the p polarization vector and gp^,r = 5t,s for /i = 0, z, 
respectively. Focusing on the decay of the F® component (p°), which alone couples to the photon, the only nonzero 
antisymmetric structure constants are and representing two independent decay channels which must be 

added incoherently. The imaginary part of the components of the self energy as computed from eq. (|4^) are then 


ImE™ (M) = ImS^, (M)(Pr + P"®) 

ImE0;(M) = ^ImE^,(M)ImE^,(M)(P°* + P°*) 

ImEtiM) = ImE^,(M)(P2^ + P^-') , (49) 


with 


ImE^^(M) 

ImE^,(M) 


327r 

327r 


(50) 


We now see that the self-energy is not transverse ^ 0) unless gs = gr- This observation enforces the equality 

of the two couplings (denoted henceforth by g), and therefore the equality of the self-energies (henceforth denoted 
by ’SiTj'Tr). The self-energies are resummed through a Schwinger-Dyson equation, yielding longitudinal and transverse 
p-propagators of the form 


DL{M,q) 
Dt{M, q) 


1 

ikf zTT,^ T Tiirp j^2 ffeET,.^,. zlmE^j.^,. 
1 

— ReE^^ — zIitlE^tt 


(51) 

(52) 


The real parts of the self energy are more involved and could be obtained from subtracted dispersion relations. We 
will assume that they do not significantly alter the pole structure. Combining eqs. (H and (1^ ) provides an 

expression for the thermal emission rate (as a function of the four-momentum g^) from the p contribution. Numerical 
integration over 3-momentum yields the rate dR/dM^ as a function of invariant mass, cf. Fig. where we also 
included a finite pion mass. The pion mass arises from a non vanishing current quark mass matrix and has been 
estimated at around m,r — 10 MeV to leading log accuracy within a Bethe-Salpeter approach for the pionic bound 
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FIG. 9. Invariant-mass spectrum from the longitudinal (dashed) and transverse (dotted) modes of propagation of the p 
meson in the CFL phase (pion mass 771-^ — 10 MeV) for pq = 350 MeV, T — 80 MeV and Pp-n-K set to the vacuum value of 6.0. 
The broad peaks of the longitudinal and transverse mode are distinct and discernible around 475 MeV. 


state energy (^. It simply introduces an extra phase space factor of (1 — into the imaginary part of 

the self-energy, which leads to a notable suppression of the rate only for invariant masses M < 50 MeV. Beyond that, 
the dependence of the (imaginary part of the) self energy causes a slight downward shift of the transverse p peak 
below its nominal mass at M = ms (this feature might not persist upon inclusion of corrections to the real part 
of Sy). A close inspection of the denominator in the longitudinal propagator reveals that it also peaks somewhat 
below M = ms in the invariant-mass spectrum. For very large M{> 1 GeV), the ratio of the rates (transverse to 
longitudinal) approaches a constant value of two, as expected. Numerical values for ms and mr follow from eqs. d^ ) 
and (III) with VMD enforced [ar = as = 2), and Ft given to leading logarithm accuracy by [|^ 

7^2 71^(1-exp(-2a;o)) ^ 

— 2 2 I 2 ’ /K ' (bd) 

TT^ Axq -b TT^ ^/2g 

The magnitude of the generalized p-meson contribution is compared to the perturbative quark rate in Fig. The 
former dominates over the latter in the invariant-mass region up to about 0.7 GeV, emphasizing the role of composite 
(i.e., generalized mesonic) degrees of freedom. Even the one-loop CFL corrections to the quark rate as calculated in 
sect, y can compete with effective theory results only for sufficiently large pairing gaps, cf. Figs. ^ and ||. 


IV. COMPARISON TO LOW-DENSITY APPROACHES 

The evaluation of the electromagnetic emissivities in the previous two sections applies to densities large enough 
to support the formation of superconducting quark matter. More specifically, we assumed the Cooper pairing to 
follow the symmetry pattern of color-flavor locking. Towards lower chemical potentials the color superconducting 
matter will undergo a (possibly first order) transition into ordinary hadronic matter^, which, nevertheless, is likely 
to occur at densities well above that of normal nuclear matter where hadronic medium effects will be substantial. 


®For realistic values of the strange quark mass, ms — 100 — 150 MeV, the CFL phase might well turn into the 2-flavor 
superconductor (2SC) before hadronic matter is recovered; the one-loop photon self energy for the 2SC state could be calculated 
along similar lines as done in sect. whereas the HLS construction of Sect. HI does not apply due to the absence of chiral 
symmetry breaking. On the other hand, if the 2SC window on the pq axis does not exist, the possibility for a continuous 
transition from CFL to hadronic matter has been raised. Such a scenario would in fact corroborate the evidence for duality 
that we are about to discuss in terms of the electromagnetic response. 
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FIG. 10. Dilepton rate (dashed line) including the peak from the generalized p meson for pq = 350 MeV, T = 80 MeV and 
A = 10 MeV. The solid line represents the free qq annihilation rate at order Og at that temperature and density. The p peak 
enhances significantly the rate from the paired quarks. 


A natural question to ask then is, how the CFL emissivities compare to results from dense hadronic phases at low 
temperatures, which should be thought of as an extrapolation from the low density side. Similar comparisons have 



FIG. 11. Gomparison of dilepton rates in the in medium hadronic and high density approaches, suggestive of a quark-hadron 
duality picture holding in dense matter. 

been discussed before for both photon and dilepton rates at large temperatures (T = 150 — 200 MeV) with and 

without the inclusion of baryon density effects. The low temperature and high density dilepton case (for 

T = 80 MeV and ns = 3.5no corresponding to a quark chemical potential of pq = 290 MeV) is exemplified in Fig. 
where the CFL rates (the combined contribution from the previous two sections) are confronted with in-medium 
hadronic calculations based on ref. |^J^, as well as the pertinent lowest order results, i.e., perturbative 0{a'l) 
qq as well as free 7r''"7r“ annihilation. The latter two are obviously very different in shape, but upon the inclusion 
of interactions, it appears intriguing that the rates from the two in principle very different ground states exhibit a 
remarkable tendency to coincide within a factor of two (somewhat more towards very low masses below 300 MeV or 
so). This agreement opens the possibility that the dilepton signal from the baryon-dense hadronic phase smoothly 
matches the description within the plasma phase, i.e., the ‘bottom-up’ and ’top-down’ extrapolations support the 
emergence of duality. 
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FIG. 12. Photon rates from CFL quark matter and a dense hadronic environment at equivalent densities. 

Another important issue in this context is the one of chiral symmetry restoration. This obviously requires an 
evaluation of the axial vector current correlation function. Although chiral symmetry in the CFL phase is broken, it 
has been found in ref. that, to leading logarithm accuracy, the masses of the generalized vector and axial vector 
mesons are degenerate. My = Ma — 2A(1 — with C = (3 — •\/ 3 ) 7 r^-\/ 6 . This is possible due to the more 

involved symmetry breaking pattern (as compared to the usual vacuum case), which, e.g., does not require a finite 
value for the two-quark condensate, but only a nonzero chiral four-quark condensate. The leading log p-ai degeneracy 
is also consistent with our earlier statement that to this accuracy, there is no coupling of the vector current to two-pion 
states. Although the simplest version of the HLS approach (as constructed in Sect III) does not incorporate explicit 
axial vector degrees of freedom, it can be accordingly extended |^] , which will not be attempted here. We furthermore 
note that model calculations |^,|| indicate certain chiral breaking effects in the CFL ground state to be substantially 
weaker than in the vacuum. For e.g., ’constituent’ quark masses turn out to be about an order of magnitude smaller 
than the pairing gaps A ~ 100 MeV. On the other hand, the pion decay constant, which also constitutes an order 
parameter of spontaneous chiral breaking, is parametrically large, {F^ = 3 Fa_^ /t^/tt^), i.e., of order O(g^fig). Thus 
the question arises how this can be reconciled with the Weinberg sum rules . One of them relates the pion decay 
constant to the integrated difference in spectral strengths for vector and axial vector correlators, 

— 4 (Imnv(qo) - Imn^(qo)) (54) 

J % 

and has been shown to hold in the medium as well A near degeneracy of Im Hy and Im 11^ implies a small 
pion decay constant. This discrepancy has been resolved in ref. |^, being associated with genuine two-loop diagrams 
involving an intermediate gluon propagator (’dumbbell diagram’). Naively, these diagrams are of higher order in 5 , but 
an infrared singularity in the gluon propagator, being regulated upon resummation by the standard Debye-/Meissner 
masses, turns them into order 0{g°) contributions. This lifts the V-A degeneracy and restores the sum rule eq. (|^). 

Finally, turning to thermal photon production. Fig. |l^ confronts rates from the CFL phase with in-medium hadronic 
results p4[| . Here, too, the emission pattern from the dense quark system is quite similar to that computed in hadronic 
or QGP scenarios. 


V. CONCLUSIONS 

In this article we have investigated the electromagnetic response of color superconducting quark matter for three 
massless flavors in the so-called CFL realization. Two approaches were pursued to calculate the e.m. current correlator. 

The first one consisted of a perturbative one-loop calculation in terms of the appropriate quark (hole) quasiparticles. 
In contrast to the free Fermi gas case, where imaginary parts arising from Landau damping are purely spacelike, the 
presence of finite quark gaps induces a non-vanishing strength (imaginary part) in the timelike region, thus opening up 
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the possibility for dilepton production. More quantitatively, taking as a baseline the free quark-antiquark annihilation 
process, quark gaps on the order of 50 MeV lead to appreciable contributions in the 3-momentum integrated dilepton 
rate up to invariant masses as large as 400 MeV. 

The second approach was based on effective field methods utilizing the notion of spontaneous chiral symmetry 
breaking in the CFL ground state. After identifying the suitable set of eight generalized Goldstone bosons, a chiral 
Lagrangian has been constructed with massive vector mesons introduced using the hidden local symmetry technique. 
Vector meson dominance arises as a natural parameter choice within this framework. The resulting dilepton rates, 
corresponding to the annihilation of generalized pions via generalized p mesons, receive contributions from both the 
spin-transverse and -longitudinal modes, the former being dominant at all invariant masses. Both modes have a 
typical mass around 500 MeV. 

Furthermore, we argued that the absence of a pion coupling to the vector current in the leading logarithm approxi¬ 
mation implies that both perturbative and effective Lagrangian contribution add to the e.m. correlator. Extrapolating 
the obtained dilepton rates down to moderate densities, we find them to be comparable (within a factor of 2 or so) 
with in-medium hadronic calculations as available in the literature. The two generic features shared in both pictures 
are (i) a strong low mass enhancement originating from soft many-body excitations, and (ii) a rather structureless 
’continuum’ beyond invariant masses of about 0.5 GeV. This extends the notion of ’quark-hadron duality’ into the 
baryon-dense regime at low temperatures, similar to what has been suggested before for the high temperature case. 
Our analysis can and should be extended to the 2SG phase, as well as the crystalline phases of dense QGD at low 
temperature |^,Q. The results should prove useful for a way to probe and discriminate the features of the QGD 
phase diagram at low temperature and moderate densities in a way that complements detailed hadronic calculations. 
Optimistically, one may even hope to use electromagnetic signals in future heavy-ion experiments for a potential 
access to the high density regime at moderate temperatures. 
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